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TECHNICAL NOTE 287 k 

ON TRAVELING WAVES IN BEAMS 
By Robert W. Leonard and Bernard Budiansky 


SUMMARY 


The basic equations of Timoshenko for the motion of vibrating non- 
uniform beams, which allow for effects of transverse shear deformation 
and rotary inertia, are presented in several forms, including one in 
which the equations are written in the directions of the characteristics. 
The propagation of discontinuities in moment and shear, as governed by 
these equations, is discussed. 

Numerical traveling-wave solutions are obtained for some elementary 
problems of finite uniform beams for which the propagation velocities 
of bending and shear discontinuities are taken to be equal. These solu- 
tions are compared with modal solutions of Timoshenko's equations and, 
in some cases, with exact closed solutions. 


INTRODUCTION 


The theoretical analysis of transient stresses in aircraft wings 
and fuselages subjected to impact loadings has generally been performed 
by means of a mode -superposition method that uses the natural modes of 
vibration predicted by the elementary engineering theory of beam bending. 
(See, for example, refs. 1 to 3.) For very sharp impact loadings, how- 
ever, this approach is known to have certain shortcomings. For sharp 
impacts of short duration, many modes are often required for a satis- 
factory degree of convergence (see, for example, ref. k) } in addition, 
the use of elementary beam theory in the calculation of the higher modes 
of vibration is inaccurate because of the neglect of, among other factors, 
the effects of transverse shear deformation and rotary inertia which 
become increasingly important for higher and higher modes (ref. 5)» 

A classically recognized alternative to the modal method of calcu- 
lating transient stresses in elastic bodies is the travel ing-wave approach, 
which seeks to trace directly the propagation of stresses through the 
body (ref. 3). Although the traveling-wave concept has been successfully 
used to treat such simple problems as longitudinal and torsional impact 
of rods, only recently have serious attempts been made to study the 
transient bending response of beams by this approach. Flugge (ref. 6) 
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was apparently the first to point out that elementary beam theory could 
not serve as an adequate basis for the traveling -wave approach since the 
elementary theory predicts that disturbances propagate with infinite 
velocityj he showed, however, that a traveling-wave theory could be 
constructed by modifying the elementary theory, as Timoshenko (ref. 7) 
did, to include first-order effects of transverse shear deformation 
and rotary inertia. On the basis of this more accurate theory, Flugge 
found that discontinuities in moment and shear travel along the beam 
with finite, and generally distinct, velocities. A similar analysis 
was carried out independently by Robinson (ref. 8 ) who, exploiting the 
hyperbolic nature of Timoshenko's equations, proposed the use of approxi- 
mate methods of solution and gave some numerical results for a particular 
example. Pfeiffer (ref. 9) also suggested the possibility of step-by- 
step solutions by the method of characteristics. In reference 10 , Uflyand 
attempted an analytical solution of Timoshenko's equations for the case 
of a simply supported beauj subjected to a sudden application of load} 
however, as was shown by Dengler and Goland (ref. 11 ), Uflyand 1 s work is 
marred by the fact that he applied boundary conditions that are incorrect 
for Timoshenko's theory. The only known example of an exact traveling- 
wave solution based on Timoshenko's theory was presented by Dengler and 
Goland for the case of an infinitely long beam subjected to a concentrated 
impulse . 

Thus, although the use of Timoshenko's theory as a basis for the 
transient -stress analysis of beams has been seriously considered, few 
problems have actually been solved. Additional basic studies of 
Timoshenko's equations and their solution, particularly for finite- 
length beams, constitute necessary prerequisites to the development of 
practical methods of dynamic -stress analysis based on the traveling -wave 
approach. In the present paper, several specific problems of transient 
loading of uniform beams of finite length are considered and their solu- 
tions by various procedures, all based on the Timoshenko theory, are 
presented. These procedures are (a) numerical step-by-step integration - 
the "method of characteristics," (b) mode superposition, and (c) exact 
closed-form solution. The examples are all for the special case of equal 
propagation velocities of shear and bending disturbances} only for this 
case have exact solutions been found in closed form. For the sake of 
completeness, the presentation of these solutions is preceded by an 
exposition of the basic equations of Timoshenko ' s theory, a derivation 
of the characteristic lines and characteristic forms of these equations, 
and a discussion of their implications concerning_propagation of 
disturbances . 


SYMBOLS 


A 


cross-sectional area 


E 


Young's modulus of elasticity 
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G 

I 

L 

M 

M 

R = 

V 

V 


shear modulus of elasticity 
cross-sectional moment of inertia 

length of finite beam (arbitrary length for infinite beam) 
internal bending moment (see fig. 1) 
dimensionless internal bending moment , ML/Elg 



vertical shear force on a cross section (see fig. 1) 
dimensionless vertical shear force, VL^/ Elg 


C 1 


c 2 


P 

q 

q 

r 

t 

v 

v 

X 


y 


y 

5 


2 ri 


propagation velocity of bending discontinuities, 

propagation velocity of shear discontinuities, 

operator used in the Laplace transformation 
intensity of distributed external loading 
dimensionless intensity of external loading, qL~’j EIp 
cross-sectional radius of gyration 
time 

velocity of deflection, y^ 
dimensionless velocity of deflection, v/C 2 
coordinate along beam 
deflection (see fig. 1) 
dimensionless deflection, y/L 

denotes a discontinuity in quantity immediately following 
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| dimensionless coordinate along team, x/L 

p density of beam material 

t dimensionless time, c^t/L 

i|r rotation of cross section about neutral axis (see fig. 1) 

a> velocity of rotation of cross section, tt 

cd dimensionless velocity of rotation, ccL/cp 

Subscripts: 

B contributing to resistance of beam to bending 

i contributing to inertia 

S contributing to resistance of beam to shearing 

x,t,i,T indicate partial derivatives with respect - to those quantities 

BASIC EQUATIONS 

Nonuniform Beams 

Timoshenko's equations .- The Timoshenko theory of~beam bending 
(ref. 7 ) constitutes a modification of elementary beam theory that' 
attempts to account for the effects of transverse shear deformation and 
rotary inertia; the basic assumption of elementary theory - that" plane 
sections remain plane - is retained, . The moment M, shear V, deflec- 
tion y, and cross-sectional rotation ijr of a nonuniform beam subject 
to a dynamic lateral loading q are governed, according to this theory, 
by the following four simultaneous partial-differential equations (see 
fig. 1): 


M +EI b * x = 0 

(la) 

v - A s G(y x - *) = 0 

(lb) 

M x - V + pli* tt = 0 

(lc) 

v x - pAiytt + 3 = 0 

(Id) 
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The first two equations constitute elastic laws relating the deformations 
to the internal loading. Equation (la) expresses the same relationship 
between moment and cross-sectional rotation as that given by elementary 
beam theory. Equation (lb) stipulates a linear relationship between the 
shear V and the shear angle y x - i|r at the neutral axisj Ag is the 
so-called "effective" shear -carrying area, different from the total 
area Ai since the true shear angle actually varies over the cross 
section. Equations (lc) and (id) prescribe rotational and translational 
equilibrium, respectively, with the term pli^tt representing the con- 
tribution of rotary inertia. 

The moment and shear may be eliminated from equations (l) to yield 
two simultaneous partial-differential equations in y and 


( ei B* x ) x + V^x “ *) ” P^tt " 0 


&S G (yx - “ P^tt + 4 = 0 


( 2 ) 


This form is convenient for finding the normal modes and frequencies of 
free vibration (q = 0) predicted by Timoshenko's theory and for carrying 
out modal analyses that make use of these modes. In this theory, a 
natural mode is described by a pair of functions [y(x) ,iKx)[] rather 
than a single function y(x), as is the case in elementary beam theory. 

For a traveling -wave analysis, however, it. is advantageous to return 
to the original system of equations (eqs. (l)) but to put them in a more 
convenient form by differentiating equations (la) and (lb) with respect 
to time. The equations become: 



(3a) 

- *55 v * - “ - 0 

(3Tj) 

+ pliC&t - V = 0 

(3c) 

- pAiV t + q = 0 

(3d) 


where the new variables, linear velocity v and angular velocity o>, 
have been introduced for y^ and ijr^, respectively. 
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Equations (3) comprise four first-order linear partial -differential 
equations in the four variables v, as, M, and V. Furthermore, equa- 
tions (3a) and (3c) contain derivatives of only M and a>, whereas 
equations (3b) and (3d) contain derivatives of only v and V. These 
facts are exploited in the next section in seeking characteristic lines, 
add characteristic forms, of these equations. 

Characteristics and the characteristic form of the -equations . - In 
equations (3a) and (3c), the variables "M and a> are differentiated 
with respect to both space and time) it would be advantageous to replace 
them by equivalent equations each involving only total derivative § (or 
differentials) in a particular direction in the space-time plane. The 
lines in the space -time plane having these particular directions - the 
characteristic lines or so-called "characteristics" - and the equivalent 
equations written in these directions are found as follows (ref. 12). 

A linear combination of equations (3a) and (3c), 



+ EI^ 


+ ^ + pl^ 


V = 0 




is formed, where the function [i is to be determined in such a way 
that the partial derivatives in equation (4). combine to give total 

derivatives f| and g in the direction of sn as yet uninovn charac- 

teristic line (i(a),t(a)] . In order— that the terms involving deriva- 
tives of M combine in the form 


M — + M+ — 
^ da % da 


3M 

da 


the function p must satisfy the following equation: 


dt 

- d- g _ At 

^ " dx “ dx 
da 

J 4_ 

where — is the required slope of the characteristic line. Similarly, 
dx 

in order that the terms involving derivatives of ca combine In the form 
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the function p must also satisfy the equation 


at 

pI i _ do_ _ at 
pEI B ax ~ dx 
an 


Since the characteristic slope must he the same in "both cases, p is 
aefinea by 


Pli 

ei b 


Thus, the two values of p an! the corre sponding characteristic 

n at 

slopes -r— are: 



at = 
ax Cl 



(5) 



>* 


( 6 ) 


where c B 
tions ( 5 ) 


/^B 

IfPli 

yieias 


Then, multiplying equation (4) by at ana using equa- 


— dM + pL dn> - V at = 0 
C 1 1 


at 



when 
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Similarly, using equation (6) gives 

- — dM + pl-t do> - V dtr-= 0 
C 1 1 


when 


dt = 



In an analogous fashion, it can be shown, from_ equations (3b) 
and (3c0j that 

i dV - pAi dv + fq + ~ dt = 0 


when 


dt = — dx 
c 2 

and 


1 / AgG 

- — dV - pAi dv + f q - 


when 



dt = 


_ 1 -- 

c 2 


dx 


/AcjG 

where Co =</ . 

2 V^i 

Then, the system of equations (3) has associated with it four- real 
characteristic directions and is thug "totally hyperbolic" (ref. 12). 

A network of characteristics can readily be constructed without prior 
knowledge of the unknowns M, V, o>, and v since their slopes depend 
only on the material and geometrical properties of the beam. For uniform 
beams, as well, as for tapered beams having uniform material properties 
and geometrically similar cross sections, cj_ and C 2 are constanty 

and the characteristics will therefore be straight; in general, however, 
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the characteristics may he curved. Figure 2 illustrates the four char- 
acteristics passing through a point P in the space- time plane with the 
characteristics designated as 



dt _ _ _1_ 
dx ~ c i 


11+ : 


dt _ _1_ 
dx Cg 


II- : 


dt = _ _1_ 
dx eg 


It is known that, by virtue of the totally hyperbolic character of the 
basic equations, the values of the unknowns M, V, cn, and v at the 
point P depend only on their initial values at t = 0 between the 
points x x and x 2 on the beam (ref. . 12 ). Furthermore, these values 

at P can, in turn, have influence only on pointB lying in the region 
above P enveloped by the 1 + and I- characteristics through P. 
Thus no signal can proceed along the beam with a velocity greater than 
(which is generally larger than C2j as illustrated in fig. 2). 


For the sake of easy reference, the four characteristic differential 
forms of the basic equations are grouped below. 


Along 

I+: 


-L_ 

dM + pl ± 

Along 

I-: 


1 

C 1 

dM - plj_ 

Along 

11+ : 

_1_ 

c 2 

dV - 

pAf_ dv + 

Along 

II- : 

_1_ 

Co 

dV + 

pAj_ dv + 


■2 


to - V dt 

= 0 

( 7 a) 

to + V dt 

= 0 

( 7 b) 

^pA 1 c^n + 

q^ dt = 0 

( 7 c) 

^pA ± Cgp - 

q^dt = 0 

(7<a) 


Related forms of these characteristic equations have been written by 
Robinson (ref. 8) and Pfeiffer (ref. 9). 
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Propagation of discontinuities .- Characteristics are lines across 
which discontinuities may exist (ref. 12) j indeed, this property is 
often taken as the basic definition of a characteristic. In the present 
problem, discontinuities (or jumps) in M and cd can therefore exist 
across the 1+ and I- characteristics, and discontinuities in V 
and v can exist across the 11+ and II- characteristics. Hence, 
a jump in M or to will propagate with the velocity cj_, whereas a 
jump in V or v must proceed with the velocity eg. It should be 
remarked that such discontinuities would appear only in the limiting 
case of a beam subjected to an instantaneous loading. The solution of 
such idealized problems, which are often instructive, requires a knowl- 
edge of the laws governing the variations in strength of these discon- 
tinuities as they propagate through the beam. These laws are determined 
below for nonuniform beams for which it is assumed that the condition 
C 1 = c 2 ^ oes not hold over any portion of the beam> in other words, 

the I and II characteristics are distinct. The special case where 
Ci = C 2 is considered subsequently when uniform beams are discussed. 

Let a and b be two points on a I- characteristic on either 
side of a particular 1+ characteristic 



If M is discontinuous across the 1+ characteristic, then M a - 

retains a finite value 8M as a and b are allowed to approach the 
1+ characteristic from either- side. Consequently, from equation (7b) 
written along the I- characteristic. 


— 5M - pit 8o> = 
ci H 1 


since dt approaches zero as a and b approach each other. Thus, 
everywhere along a 1+ characteristic, jumps 8M and 8 cd across this 
characteristic are related by 


8M = cj_plj_ Scd 


( 8 ) 
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Similarly, the jumps across the other characteristics can be readily 


shown 

to satisfy 



Along 

I-: 

SM = -c-^pl^ Scd 

( 9 ) 

Along 

11+ : 

SV = -c 2 pA^ Sv 

(10) 

Along 

II- : 

S V = CgpA^ Sv 

(n) 


A jump in M is thus always accompanied by a definite jump in to; 
similarly, jumps in V and v are always coupled together. 

The variations in the magnitude of a discontinuity as one proceeds 
along a characteristic may be determined in the following manner. Equa- 
tion (7a) is written for the upper side and the lower side of the 
1+ characteristic j then, since V is continuous across a 1+ char- 
acteristic, the difference of the two equations yields 


d(BM) + c-jpli d(5o>) = 0 

along I+. Eliminating 8o> by using equation (8) gives 

d(W) * c lP I ± - 0 


By carrying out the indicated differentiation in the second term and 
dividing by 2SM, the following result is obtained: 


d(5M) = l. tj(ciPlj) 
SM 2 c-jpl^ 


Solution of this equation gives 


(SM) 2 = (5M)i 


(ciPlj) 2 
(°l pI i)l 


( 12 ) 


as the relationship between the magnitude of the jumps in M at two 
points 1 and 2 on the 1+ characteristic. 
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It can be shown that the identical relationship holds between 
jumps in M at two points on a I- characteristic. Similarly, it can 
be found that, on 11+ and II- characteristics. 


(sv) 2 



(13) 


for any two points 1 and 2 on a given characteristic. The corre- 
sponding variations of the jumps Sen and &v are, of course, readily 
determined from equations (8) to (11). 


Uniform Beams 


Nondimensional form of the equations . - The examples to be presented 
in this report are all concerned with beams having uniform cross-sectional 
size and shape and uniform material properties throughout their length. 

For such beams, it is convenient to express Timoshenko's equations in 
nondimensional form~ Thus, equations (3) may be written In terms of 


M = 


ML 


TT VL" 


El- 


B 


V = 


EV 


— chL 
CD = , 


and 


v =» 


V 

c 2 


as 


CD| + My = 0 


c 2 _ 1 — _ 

_ Vg - g V T - CD - 0 


(14a) 

(14b) 


M^ + ac^ - V = 0 


h - R v r + a = 0 


(14c) 

(l4d) 


where 


t = * 

5 TV 


T = 


Cjt 


4 *I B ' 




The quantity L 


refers to the beam length for all beams except those of infinite length, 
in which case any convenient arbitrary length may be chosen for L. 


,The characteristics of Timoshenko's equations for a uniform beam 
are defined in the £,t plane by the families of straight lines 
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I+: 

*1=1 
d 4 

( 15 a) 

I-: 

*1 = _1 
di 

(15b) 

11+ : 

dT C 1 

d i “ C2 

( 15 c) 

II-s 

dr _ C 1 

d| " c 2 

(15d) 


The nondimensionalized characteristic forms of. the basic equations become 


Along 

I+: 

dM + 

to 

- V dT = 0 

, (l6a) 

Along 

I-: 

<3M - 

dZ5 

+ V dT = 0 

( 16 b) 

Along 

11+ : 

dV - R dv 

+ 1 

f c 2 \ 

EZ5 + — — q JdT =0 

^ c 1 / 

( 16 c) 

Along 

II-: 

dV + R dv 

+ 1 

fa ~ If ^) dT " 0 

(l6d) 


In addition to the restriction to uniform beams, for which c-j_ 
and C 2 are constants, the examples presented herein are further limited 
to those beams for which the propagation velocities cp and c 2 are 

equal. This assumption has been made because the simplifications that 
result not only permit the ready attainment of numerical solutions but 
also, in particular cases, permit the attainment of exact closed solu- 
tions for comparison. Since, for this very special case, the charac- 
teristics II coincide with the characteristics I, equations (l6) may 
now be written, for "5=0, as 


Along I+: 


fdM + dZ5-VdT=0 

dV - h\ 2 dv + dT = 0 


(17a) 


Along I-: 


f dM-doS + VdT = 0 

dV + hX 2 dv + dT = 0 


(1Tb) 



where 
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Propagation of discontinuities when cq = C 2 '.- Equations (8) to (13), 

which describe the behavior, of discontinuities in a nonuniform beam, also 
describe, as a special case, the behavior of discontinuities in a uniform 
beam for which cq ^ 02 - They show that such discontinuities propagate 
with constant magnitude. 

However, when the beam has properties such that cq = C 2 * these 
equations are no longer valid. Equations which are valid in this case 
may be derived in precisely the same way by using equations ( 17 ) instead 
of equations (7). Discontinuities in such beams can be shown to be 
related by 


Along I+: 



(18a) 


Along I-: 


{ 6M = -6cd 
SV = hX 2 Sv 



(18b) 


as they propagate, and they can be shown to vary in magnitude according 
to the equations 


Along 1+ and I-: <! 


i(®) - § SV - 0 


■^(SV) + 2\ 2 5od = 0 


(19) 


Equations (19) may he solved simultaneously to obtain 


Along 1+ and I-: 


f Scd = A cos(Xr - B) 
\_6V = -2XA sin(XT - B) 


( 20 ) 


where A and B are arbitrary constants which must be evaluated by 
using known values of Sco and SV at some point’ on the characteristic. 
The variations in 5M and 8v can then be readily found by using 
equations (18) . 
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Thus, for the case c-j_ = c^ } discontinuities in a uniform beam do 
not propagate unchanged but vary in magnitude sinusoidally as they 
progress through the beam. 


Limitations of the Theory 

It may be well to insert a word of caution about the applicability 
of Timoshenko's theory. The investigations of Prescott (ref. 13) and 
Cooper (ref. 14) have shown that, when the response of a beam includes 
components whose wave length is small compared to the depth of the beam, 
the assumption of Timoshenko's theory that plane sections remain plane 
after bending is, as might be expected, too restrictive. Since applied 
disturbances which could give rise to discontinuities would obviously 
excite even the shortest wave length in the spectrum of the response, 
the results obtained by application of Timoshenko's theory to such 
hypothetical problems cannot, in themselves, have practical significance. 
However, carrying out solutions involving discontinuities is a useful means 
of testing methods of solution of the Timoshenko equations with a view to 
applying these methods to problems in which discontinuities do not exist. 
Furthermore, the admittedly inaccurate response to an infinitely abrupt 
disturbance may be used to obtain the correct response to disturbances 
of a more practical nature through the application of Duhamel's super- 
position integral. 


SPECIFIC EXAMPLES - FINITE UNIFORM BEAMS WITH c ± = c 2 


Three specific examples are considered; they are: a cantilever 

beam given a step velocity at the root, a simply supported beam subjected 
to a step moment at one end, and a simply supported beam subjected to a 
ramp-platform moment at one end. 


Methods of Solution 

In the examples, the results of calculation by the following three 
methods are compared: (a) numerical step-by-step integration along the 

characteristics, (b) normal -mode superposition, and (c) exact closed- 
form solution. The first two procedures are approximate in character, 
but they could conceivably be generalized sufficiently to be applied to 
practical structures; the last procedure, although exact, would rarely 
be useful in practice and Is introduced herein chiefly as a check on the 
accuracy of the first two. 

The detailed descriptions and applications of the three methods are 
contained in appendixes A, B, and C. In brief, the numerical procedure 
exploits a grid of characteristic lines as shown in figure 3(a) . For 
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the case cp = c 2 that is under consideration, this grid consists of 

two families of" lines in the | ,t plane with slopes 1 and -1. A 
recurrence formula is developed in appendix A that gives the values 
of ci5 and V at station 1 of a typical interior mesh (see fig. 3(h)) 
in terms of the values of 55 and V at stations 2, 3 > and 4. Repeated 
application of this formula, together with the use of special formulas 
for the half -meshes at either end of the beam and the knowledge of-the 
magnitudes of jumps in 55 and V across characteristics where they 
occur, permits and V to_be calculated throughout- the |,t plane. 
Subsequent determination of M and v is achieved by means of .simple 
addition formulas utilizing these calculated values of 53 and V. 

Although the solutions derived in appendix B have actually been 
obtained by Laplace transform techniques, they have been termed "modal 
solutions" because they axe exactly those which would result from an 
application of the usual mode -superposition process. The exact closed 
solutions in appendix C have also been obtained through the use of 
Laplace transforms. 


Cantilever Beam Given a Step Velocity at the Root 

The first example to be considered, the response of a uniform 
cantilever beam given a step velocity v = 1 at the root, is the 
equivalent of the most severe idealized landing problem, the instan- 
taneous arrest of the root of a moving cantilever beam. Computed results 
for the shear and moment at' the root -of such a beam having a ratio of 
length to radius of gyration of 10 (A, = 5) and properties such that 
C 1 = c 2 8X6 presented in figures 4(a) and 4(b), respectively. Two 

separate curves obtained by the numerical procedure are shown - one from 
a grid that divides the beam into 10 segments and the other from a 
20-segment solution. The modal solution includes the contributions of 
the first eight modes. Results given by an exact closed solution are 
shown for both the shear and moment at the root up to the time- r = 2. 
These exact results are actually those for an infinitely long beam, since 
the influence of the free end is not felt at - the root until T = 2. 

After r = 2, the influence of the free end is felt' and, in this case, 
an exact solution is not feasible. To illustrate the time range covered 
in the plots, the point corresponding to half the period of the first 
mode of vibration of the beam is indicated on the time scale of each 
plot. 


In figure 4(a), the shear discontinuities evident in the numerical 
solutions occur each time the discontinuous wave front returns to the 
root after being reflected at the free end. The beam is seen to react 
violently to each of these boosts by the wave front, with more and more 
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osclUations occurring after each succeeding jump. The frequency of 
these oscillations tends to increase with each succeeding boost until 
limited by the finite time interval. In these regions of violent 
oscillation the accuracy of the numerical solutions is obviously ques- 
tionable; indeed the question arises as to whether these oscillations 
are really predicted by the theory or are merely the result of some 
instability in the numerical process. This question is resolved in the 
next section in which the simply supported beam is considered. At any 
rate, away from the regions of violent oscillations, the comparisons 
with the modal solution are favorable, and, for t < 2, the fine-grid 
numerical solution almost coincides with the exact closed solution valid 
in this region. 

The comparisons between the numerical and modal solutions are very 
good in figure 1(b) where the time history of the moment at the root is 
plotted. Again, the fine -grid numerical solution nearly coincides with 
the exact closed solution for T < 2. 


Simply Supported Beam With an Applied End Moment 

Step moment .- A simply supported uniform beam with a ratio of 
length to radius of gyration of 10 (X = 5) and properties such that 
C]_ = cp is subjected to a step moment M = 1 at the end £ = 0. Com- 
puted results for the shear at 1=0 are presented in figure 5(a) and 
the time history of moment at the center of the beam is presented in 
figure 5(l0« The point corresponding to the full period of the first 
mode Is marked on the time scale of each plot. The numerical curves 
for both shear and moment were obtained by a 20-segment solution. The 
modal curves were obtained by adding dynamic corrections to the static 
solutions, the dynamic corrections being calculated with six modes 
for both the shear and the moment. 

This example affords an answer to the question raised In the 
preceding section with regard to the stability of the numerical procedure 
following the passage of a discontinuous wave front. The fact Is that 
the violent oscillations that occur after the discontinuity actually 
appear in the exact solution (fig. 5(a)) and are hence inherently present 
in the theory. 

In figure 5(a) the scale of dimensionless vertical shear happens 
to be precisely the dynamic overshoot factor; it is of interest to note 
that values at least 15 times the static shear are predicted when the 
moment Is applied suddenly. 

For shear (fig. 5(a)), the Inaccuracies in the numerical solution 
just after the discontinuities are evident; however, the numerical 
results approximate the exact solution very well elsewhere. A similar 
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observation may be made for the modal solution, which, as would be 
expected,, ignores the discontinuities and violent' oscillations caused 
by the wave front. The numerical and modal curves for the moment in 
figure 5(b) follow the same pattern, agreeing well everywhere except 
in the regions immediately following the discontinuities. 

Ramp -platform moment . - Computed results for the shear at the end 
5=0 and the moment at the center of the same simply supported uniform 
beam subjected to the applied ramp-platform momentr 


M(0,t) = t (0 ^ t % 1) 

M(0,t) =1 (t > 1) 


are presented in figures 6(a) and 6(b). Again, the period of the first 
mode of -vibration is marked on the time scales. 

The numerical curves were recomputed with a 20-segment grid for 
the new forcing function} however, the results for the modal and exact 
solutions were obtained by means of a superposition of the preceding 
results. This superposition was carried out analytically for the com- 
plete modal solution and for the exact solution in the range t < 2. 

In the range r > 2, it was necessary to carry out the superposition 
for the exact solution numerically. 

In figure 6 the time to peak value of the applied moment is 
seen to be approximately one-seventh the period of the first mode; 
predicted shear values approximately three and one-half times the 
static response occur. 

With the removal of the discontinuity, there remain no high- 
frequency oscillations which the numerical solution might be unable 
to represent. In fact, all three solutions for the shear and both 
moment solutions are seen to be in excellent agreement. 


CONCLUDING REMARKS 


Timoshenko's equations for the motion of vibrating nonuniform beams 
may be written in a characteristic form which appears to be well- suited 
to solution by numerical methods. In the examples presented in this 
report, all of which are for uniform beams with equal propagation veloc- 
ities of bending and shear discontinuities (cp = cq) > the solutions by the 
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numerical and modal methods generally agree well with each other as well 
as with exact closed solutions where these have "been obtained. However, 
the modal method, of course, fails entirely to reproduce discontinuities 
in shear or moment, and the numerical procedure, although it yields these 
discontinuities, does so with decreasing accuracy as more and more 
reflections of the wave front occur. In the more practical situations 
where discontinuities do not exist, these difficulties will, of course, 
not arise. 

The results of the examples carried out by the numerical traveling- 
wave procedure encourage the viewpoint that such traveling- wave analyses 
may eventu all y be of practical usefulness. This kind of procedure is 
inherently simple and straightforward and has the advantage that the 
bulk of the labor involved is routine computation rather than math- 
ematical analysis. It should be emphasized, however, that numerical 
solutions of Timoshenko's equations have been demonstrated only for 
uniform beams in which the propagation velocities C]_ and eg a re 
equal} numerical procedures for the general case where they are unequal 
remain to be developed and tested. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 28, 1952. 
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APPENDIX A 

NUMERICAL SOLUTIONS FOR UNIFORM BEAMS WITH c ± = c 2 
Matrix Formulation 


Let the differential equations (17) tie replaced "by the finite- 
difference equations 


Along 1+ : 


Along I-: 


AM + Aui - V At = 0 

(Ala) 

AV - 4X 2 Av + 4X% At = 0 

(Alb) 

AM-Aj3 + VAt = 0 

(Ale) 

AV + 4X 2 Av + 4X 2 S At = 0 

(Aid) 


and consider a closely spaced network of 1+ and I- characteristics 
in the space-time plane as shown in figure 3(a). Let the corners of a 
typical interior mesh of this grid be designated_as shown in figure 3(t). 
A step-by-step integration formula for ca and V may now be derived in 
the following manner. 

Equations (Al) may be written along the upper sides of the typical 
mesh to obtain 


— M 2 + oo^ — ciig 

- + v 2 ) = 0 

(A2) 

M 1 ~ M 3 - o>l + 0 D 3 

+ ^(v x +.v 3 ) = 0 

(A3) 

_ v 2 - 4x 2 (v 1 - v 2 ; 

) + 2X 2 At^53l + 53^ = 0 

(A4) 

- V 3 + 4X 2 (v x - v 3 ; 

) + 2X 2 At^5>[_ +53 3 ) = 0 

(A5) 


where 53 and V have been assumed to vary linearly in the small 
intervals between the corners. Obviously equations (A2)_to (A5) may __ 
be solved simultaneously to obtain the four quantities V, 53, v, and M 
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at point 1 in terms of their values at points 2 and 3» However, it is 
noted that, if V, 35, v, and M at points 2 and 3 'were determined by 
a s imil ar process for the preceding meshes, they already satisfy the 
equations 

Mq - M 4 - 352 + 35^ + ^(? 2 + V 4 ) = 0 (A 6 ) 

M 3 - + S 3 - 554 - £g(v 3 + V 4 ) = 0 (A7) 

V 2 - 74 + 4l 2 (v 2 - V 4 ) + 21 2 At(35 2 + 354 ) =0 (A 8 ) 

V 3 - V 4 - ^ 2 (^3 - V 4 ) + 21 2 At^3 > 3 + 5 S 4 ) = 0 (A9) 

Equations (A 6 ) and (A7) may be added to equations (A2) and (A3), 
respectively, to obtain 

M-l - M4 + 35 ^ - 2S 2 + 554 - - V4) = 0 

and 

M-l - M 4 - 3>l + 235 3 - 534 + ^(v-l - V 4 ) = 0 

which may, in turn, be subtracted to obtain tbs single equation 

3^ - 33 2 - oi 3 + Q 4 - - V 4 ) =0 (A10) 

Similarly, equations (A 8 ) and (A9) may be subtracted from equations (A4) 
and (A5), respectively, to obtain 

Vi - 2V 2 + V 4 - 4l 2 (vj_ - V 4 ) + 2X 2 At^55 1 - CD 4 ) = 0 

and 

Vt_ - 2V 3 + V 4 + 4X 2 (v x - V 4 ) + 2X 2 At(33 1 - © 4 ) = 0 
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and these may he added to derive 


V x - V 2 - V 3 + V 4 + 2X. 2 Ar(c5 1 -55^) =0 (All) 

Equations (A10) and (All) may now be solved simultaneously to obtain 
the step-by-step integration formula 


“1 

Vl 


1 . 

(X At) 2 + 1 




(A12) 


where 


and 



At 

2 

1 



|(X At) 2 - 1 
4x 2 At 


-At 

(X At) 2 - 1 


Thus E and V may be determined at every interior mesh point 
in the space-time plane by the repeated application of formula (A12) to 
each mesh as it is encountered. The half -meshes which occur at' the 
vertical left and right boundaries of the plane (fig. 3 (a)) require 
special formulas incorporating the known boundary data. These formulas 
may be derived from equations (Al) by a procedure similar to that used 
in obtaining equation (A12). Boundary formulas for some specific 
examples are presented in subsequent sections of this appendix. 


Besides boundary data, initial data must be provided in order that 
a step-by-step solution may be begun. In all the examples considered, 
the beam is Initially at rest and has a disturbance applied at the 
point |=0 beginning at time t = 0. The region t < % (without grid 
lines_in fig. 3(a)) therefore is one of zero stress and motion and 35 
and V are given along the line r = | by the known conditions at the 
wave front. It is with these values that each numerical solution is 
begun. 
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Discontinuities in 55 and V offer no special difficulties since 
they propagate so as to he always located on I characteristics which, 
of course, can he made part of the basic grid. Thus they are simply 
added as they are encountered. 

Once 55 and V have heen determined at the grid intersections hy 
the step-by-step process, M and v may he found hy repeated applica- 
tion of equation (A2) or (A3) and equation (A4) or (A5) along the proper 
grid_lines from boundaries where M and v are known. For example, 
if M is known along the left boundary, it may he found at a point g 



hy applying equation (A2) successively to the__intervals ah, he, 
cd, . . . fg. The resulting expression for M g becomes 

Mg = M a + Eg. - 55g + At(| V a + + . . . + V f + | V g ) (A13) 


This procedure is seen to correspond to integration of the first of 
differential equations ( 17 a) by means of the trapezoidal rule. 

Specific Problems 

Cantilever beam given a step velocity at the root . - If the 
root 1=0 of a uniform cantilever beam is given a step velocity v = 1 
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at time r = 0 , the boundary conditions may be written 


v(0,t) = 1 
E(0,t) = 0 
M(l,r) = 0 
V(l,r) = 0 






J 


(A14) 


Typical left- and right -boundary half -me shea are 



Application of equations (Alb) and (Aid) to the sides of the left - 
boundary half -mesh and proper combination of the resulting equations 
to eliminate v produces 


Vi = 2 V 3 - V 4 (*A15) 

Similarly, application of -equations (Ala) and (Ale) to a right-boundary 
half -mesh produces 


= SEBg - coi. 


(A16) 


so that a complete set of integration formulas for determining £3 and V 
for this problem is now available. In addition, if- the unused equations 
(eqs. (Ala) and (Ale) for the left-boundary half-mesh and eqs. (Alb) 
and (Aid) for the right) are applied and combined to eliminate K and v 
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at points 3 and 2 , the recurrence relations 


M x = M 4 - - ~{V 1 - V 4 ) 


(AIT) 


for left -boundary half -meshes and 

7 1 = ^ " ^2 ^2 + " “4) ( Al8 ) 

for right -boundary half -meshes are obtained. These equations may be 
used to compute M at the root and v at the tip after 53 and V have 
been determined everywhere. 

Since the applied disturbance is initially discontinuous, 53 and V 
will be discontinuous along the line t = | (see_fig. 3(a)). In order 
to begin the step-by-step solution for 53 and V, these discontinuities 
must be determined in advance. Furthermore, since the wave front is 
reflected back into the beam at either end, 55 and V will also be 
discontinuous along the lines t = 2 - t = 2 , and 655 

and 8 V must be determined along each of these lines before the step- 
by-step solution can be extended beyond it. The discontinuities at the 
wave front are determined as follows. 

From the boundary conditions (Al4) , it is seen that 655(0,0) = 0 

and 8v(0,0) = 1. Thus, from equation (l 8 a), 
equations ( 20 ) become, for the line t = |, 

8o3(£,t) = -21 sin 1 | 

8 V( |, |) = -4l 2 cos X! 


5V(0,0) = -4l 2 , and 


(A19) 


At |=1, the discontinuities across t = | are 605 ( 1 , 1) = 5M(l,l) = 
-21 sin_l and 6V(l,l) = -4l 2 cos 1, so that, if the boundary condi- 
tions M(l, t) = v(l, t) = 0 are to be satisfied, jumps must occur across 
t = 2 - | at the point (1*1) with the magnitudes £53(1,1) = -6M(l,l) = 
-21 sin 1 and 5V(l,l) = 4l 2 cos 1. In view of these initial condi- 
tions at t = 1, equations (20) become, for the line r = 2 - |, 


So>(i,2-!) = -21 sin l£ 
5V( 1,2-6) = 4l 2 cos 1| 


(A20) 
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Initial jump values £xo(0,2) = 0 and 5V(0,2) = 4X 2 for the line 
t = 2 + 6 may be found from equation (A20) by satisfying the boundary 
conditions 55(0 , t) = 0 and v(0,t) = 1. Prom these initial values, 
the discontinuities across t = 2 + £ are found to be the negative of 
those across t = or 


Sm( |,2+|) = 2X sin X§ 
SV( I,2+|) = 4X 2 cos X| 


Then, it must be true that 


8a>( l,4-|) = 2X sin X£ 
6V(l,4-i) = -4x 2 cos XI 


(A21) 


(A22) 


and so forth, with the values on each succeeding line t = n t £ repeating 
the values on the line t = (n - 4) ± £. The variations in the magnitudes 
of S55 and SV as the wave front propagates back and forth through the 
beam are thus clearly defined, and, since t < | is a region of zero 
stress and motion, equations (A19) define the values of <jo and V along 
the line t = 

With 55 and V known along t = g, the solution for 55 and V 
may be begun by applying formula (A15) to the half -mesh in the lower 
corner of triangular region (l) (fig. 3(a)). It is continued step by 
step throughout the triangle, with formula (A12) being used for all 
interior meshes. When the first triangle is complete, the known jumps 
along the line t = 2 - £ may be added to the values computed for the 
under side of this line, and the solution may then be carried out in 

triangle (g) beginning again, at the lower corner this time with 
formula (Al6) . In this way the solution may be carried through as many 
triangular regions as desired. 

Two sets of computations, with time intervals _At = 0.1 and 
At = 0.05, have been made for a cantilever beam for_ which X = 5- These 
time intervals correspond to grids dividing the beam intor 10 and 
20 segments, respectively. (The computations for the 20-segment solution 
were made on the Bell Telephone Laboratories X-66744 relay computer at 
the Langley Laboratory. ) The computed time histories of shear and moment 
at the root have been plotted in figures 4(a) and 4(b), respectively, up 
to tj= 10. In order to obtain these time histories, the computations 
for and V had to be carried through the firstnine triangular 
regions . 



NACA TH 2874 


27 


Simply supported beam with an applied end moment . - If the end | - 0 
of a uniform simply supported "beam for which C]_ = is subjected to 
an applied bending moment Mq(t) beginning at time t = 0, the boundary 
conditions are 


M(0,t) = Mq(t) 
v(0,t) = 0 
M(1,t) = 0 
v(1,t) = 0 




(A23) 


Typical left- and right -boundary half -meshes for this beam are 



Equations (Ale), (Ala), (Aid), and (Alb) may be written along the 
sides of the left half -mesh to give, respectively. 


(Mq)i - M 3 - + cd 3 + ^(v-l + V 3 ) = 0 

% " (^0)4 + ®3 - ®4 " ^ 5^3 + V4) = 0 
V x - V 3 - 4x 2 v 3 + 2X 2 AT^c5i + oi 3 ) =0 
V 3 - V 4 - 4 X 2 v 3 + 2 cd ^ + S 4 ) = 0 



28 


NACA TN 2874 


Adding the first two equations and subtracting the last two equations 
produces 

(Mo)! ~ (Mq ) 4 - 5^ + 2 E 3 - E k + ^{Vi - V 4 ) =0 

and 


1 ± - 2V 3 + Vi,_+ 2X 2 At^o5 1 -55^) = 0 


which may be solved simultaneously for o>]_ and V]_ to obtain 
n- -r 


“1 


(X At ) 2 + 1 


A23 


y- -> 


205- 


2V- 


a 4 


014: 


1 - 


-21 2 At 


(A24) 


where and have "been defined in conjunction with equa- 

tion (A12). 

A similar process produces, for- the right -boundary half-mesh, the 
formula 



1 

“ — 


2&2 



■> 

CD4 

Vl 

— H 

(X At ) 2 +-1 

A23 


2V 2 

+ 

A4 

► 

^4 

J 


(A25) 


Now consider the case in which the applied moment is a step moment; 
that is, Mq( t) = 1. Equation (A2b) reduces to 


0>1 


(X At ) 2 + 1 


^23 

vU 


2d>z 


2V-. 


A4 


C04 

^4 




(A26) 


for left-boundary half -meshes. A similarity between the boundary 
formulas (A25) and (A2 6) and the interior formula (A12) is apparent; 
indeed, formula (A12) may be used everywhere if the sums of— the values 
at points 2 and 3 are simply replaced by twice the values at the inboard 
point of each boundary mesh. 
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For this problem, as in the preceding one, the magnitudes of the 
discontinuities £co and BV at the wave front must be determined in 
advance for every point on the grid lines defining its position. From 
the boundary conditions and equations (l 8 a) it is apparent that 
So5(0,0) = SM( 0,0) = 1 and SV(0,0) = 0. Thus, for the line r = g, 
equations ( 20 ) become 


Boi = cos Xt 
BV = -2X sin Xt 


(A27) 


Furthermore, it is found that, when the wave front is reflected at 
either of the simply supported ends, the signs as well as the magnitudes 
of the discontinuities remain unchanged, in contrast to the behavior of 
a discontinuous wave front reflected at a free or fixed end as in the 
preceding case. Thus equation (A27) must be valid for the entire zigzag 
path defining the position of the wave front. 

With the discontinuities known, the step-by-step solution can be 
begun as before. Again, it is convenient to complete each triangular 
region before proceeding to the next. 


Computations for this problem were made on the Bell computer, with 
a 20 -segment grid, for a beam with X = 5 • The quantities 35 and V 
were found in the first eight triangular regionsj in addition, the moment 
at the center was computed. The time histories of shear at the end | = 0 

and of moment at the center | = — have been plotted in figures 5 ( a ) 
and 5 (b), respectively. 

Now consider the case in which a ramp-platform moment is applied to 
the beam} that is, Mq(t) is defined by 


Mq(t) = r 

Mq(t) = 1 



(A28) 


Equation (A24) reduces to 




r 

— - 

1 


A 23 

(X At ) 2 + 1 






A k 

v h 


+ 2 An 


-2X 2 At 


(A29) 
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for the left-boundary half -meshes for which ^ 1 and reduces to 
equation (A26) for all the rest. Equations (A12) and (A25) are, of 
course, still valid for interior meshes and right -boundary half -meshes, 
respectively. 

Since 5oJ(0,0) = SV(0.,0) = 0 for this case, the solution begins 
with initial values o> = V = 0 along the line t = £ and there are no 
discontinuities to - be added. Thus the solution may progress in any 
convenient manner without regard for the position of the wave front. 

The same quantities were computed for this problem that were obtained 
in the preceding case. As before, it was assumed that X = 5 , and the cal- 
culations were performed on the Bell computer with a 20-segment grid. 

The resulting time histories of shear at £ = 0 and moment at £ = ^ 
are presented in figures 6 (a) and 6 (b), respectively. 
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APPENDIX B 

MODAL SOLUTIONS FOE UNIFORM BEAMS 


The solutions carried out in this and the following appendix make 
considerable use of Laplace transform techniques; all the Laplace 
transforms used, most of which were taken directly from references 15 
and 1 6 , are given in table I for the sake of easy reference. 


Equations 


The equations of motion (2) may be written for a uniform beam 
vibrating in the absence of external distributed loading in nondimen- 
sionhl form as 



+gg + R (y^ - *) - * TT = 0 


r ( B1 ) 


where y = ^. The dimensionless bending moment M and vertical shear V 
are given by (see eqs. (l)) 



(B2) 


V = R(y g - i) 


(B3) 


If the beam is initially at 
j 0 ) = 0 ), equations (Bl) may 
table I has been used) 


rest (y(£,0) = y T (|,0) = i|f(|,0) = 
be transformed to (where transform 6 in 


( Y t - *) fi 



Y (B4) 


+ *(*£ - $) - = 0 
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in terms of the Laplace transforms Y(g,p) = / e 2 y(g,T)dT and 

Jo 

p 00 

■f(|,p) = / e -pT i|f( ?, T )dT. The solutions of equations (b4) have the 

JO 


form 


mi 


Y(|,p) = A(p)e 
$(i,p) = B(p)e m ^ 


(B5) 


Substituting these relations into equations (b4) leads to a biquadratic 
equation in m 


4 2 

m - p 


1 + 


m 2 + 


(sjjW + B > - 0 


(B6) 


and gives the following relationship between B and A: 

.2 


m 


B = 




c 2/ 


m 


(B7) 


Let the solutions of equation (B6) be written m = iim^ and 



inig * 


Then 


and the general solution of equations (b 4) is 


Y(i,p) = cos + Cg sin m-j_| + cosh mgi + sinh m^g (B8) 
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m,v) 





sinh 



(B9) 


where Cg, C and are functions of p which must be chosen 

so that the boundary conditions are satisfied. 


Specific Problems 

Cantilever beam given a step velocity at the root.- If the 

_ C 1 

root 1=0 of a cantilever beam is given a step velocity v = — y T = 1 
at time t = 0 } the boundary conditions may be written 


y (°< T > T 

♦(0,T) = 0 

V|(1,t) = 0 

t(l,r) = yj(l,T) 


(BIO) 


These boundary conditions transform to (transform 8 , table I) 


T(0,p) 


f2 _1_ 

C 1 p 2 


$(0,p) = 0 

f|(l,p) = 0 


f(l,p) = T|(l,p) 




> 


(Bll) 


The constants 811(1 c 4 Bxe determined "by substituting 

equations (B8) and (B9) into equations (Bll)j the resulting expressions 
for Y(|,p) and $(i,p) may be written in the forms 
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1 p 2 D(p) 


(B12) 


and 


f( 6 ,p) = 


U|(i,p) 

c i p 2 D(p) 


(B13) 


where 


u y(I>p) (l - — sin !!■[_ sinh mg + ^ cos m2 cosh mgjcos n^l + 

+ m^ sin m l slnh “2 + Q cos m x cosh m^j cosh mg? + 
^ sin m.2 cosh ^ + ~ cos m2 sinh a^j sin m x 6 - 
^ siD m l cosh “2 + Q cos m x sinh mgjsinh mg| 


u f (^P) - 


- 2 - (i)V 


mg 


fmg 


sin m2 cosh mg + Q cos m2 sinh mgj^cos m^? - 


cosh 


“ 2 ^) " Q mj ^ 1 " mg 8in ra l sinh + ^ cos m i cosh mg^sin m2? + 


^ + m x 8ln m l slnh “2 + Q cos m x cosh mg^sinh mg| 


D(p) = 2 + 


(q + i)cos m2 cosh mg + ( 2 s - ^Bin m x sinh mg 


Q = 


+ (t^J P 2 


mg £ 
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The inverse transforms y(|,r) and ^(i/r) are determined by 
substitution of Y(|,p) and \jf(|,p) into the complex inversion integral 

(ref. 15). The singularities of the resulting integrands eP T Y(|,p) 


and eP T $(|,p) must therefore be examined. Although the functions m^ 
and mg are , in themselves, multiple -valued functions of p with branch 

2l/R 


points at p = 0, ii^H, and ± 


c 2 


-, it follows from a considera- 


1 - (® 2 

tion of the fundamental theorem of the uniqueness of solutions of ordinary 
linear differential equations that Y(£,p) and $(|,p) must neverthe- 


less be single-valued. The integrands e^ T Y(|,p) and e^ rr $(|,p) thus 
have no singularities other than poles. The inverse transforms y(l> T ) 
and t(l,T) will therefore be taken as the sum of the residues at the 


poles of eP T Y(£,p) and e^ T $(|,p), respectively (see ref. 15 ). 


It can be shown that no singularities occur in the numerator 
functions Uy(£,p) and U-^(£,p)j therefore, all the poles of Y(|,p) 
and $(i,p) must be introduced by the zeros of the denominator p 2 D(p) . 
Consider first the equation 


D(p) 


2 + 



”1 


cosh 


1112 




sin m^ sinh mg = 0 


(B14) 


This equation has an infinite number of roots p = ±p n , • n = 1, 2, . . . <», 
all of which are, in general, simple poles of both Y(|,p) and \jr(|,p). 

In addition, p = 0 can be shown to be a double pole of Y(|,p) but not 
a pole of $ (£,p) . 


The sum of the residues of eP r Y(|,p) at the poles p = 0 and 

p = ±P n , n = 1, 2, 3; • • • ds 




T + 


n=I 


1 

%(^Pn) c p n x 

Uy(l^-P n ) -PnT 

Pn 

D 1 (pn) 

+ D'(-Pn) ‘ e J 
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and, similarly, the sum of the residues of e pr f(£,p) at- p = +p n , 
n = 2 , 3, . . . oo, "may he written 


fU/r) 



where 


D'(P) 



However, it is seen that U Y (|,p) ='%(|,-p), U f (|,p) = U $ (|,-p), and 

D'(p) = -D'(-p) so that these equations reduce to 
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and 


yU, T ) 



+ 




Uyd^Pn) 

Pn^’CPn) 


sinh p n T 

_ 


c p u $(t*Pn) 

t(i, T ) = 2 — ) — sinh p n r 

5=1 p n D ’(Pn) 


(B15) 


(Bl6) 


If , in equation (Bl ^) , p 


is replaced by ik, the result may be written 


2 + 



2 


cosh a cos 


3 




sinh a sin p = 0 


where 


(B17) 



and k = ft being the circular frequency of vibration. Equa- 

C 1 

tion (B17) is the frequency equation for a uniform cantilever beam (see 
ref. 5) and its solutions k = k^ n = 1, 2., 3, . . .a, are the non- 
dimensional natural vibration frequencies of such a beam. Equations (B15) 
and (Bl6) may be written in terms of these natural frequencies simply by 
replacing p n by ikjj. The results are 
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T + 2 Xy n (|)sin k n T 
n=l 


(Bl8) 


and 


t(l, T ) = 2 — X^ n (^)sin k n T 


(B19) 


n=l 


where 


X 7n (5) = |^1 + ^ sin sinh csq +-^ cos 3 n cosh cc^cosh o^t 

( 1 - — sin (3 n sinh Oq + 0n cos p n cosh c^jcos |3 n £ - 
cosh sin J3 n + sinh cxq cob p n ^sinh an I + 

sinh an cos J3 n + 0n cos ^ On s ^ n Pn^ s ^- n ^ s ] 


- (%)V 

X *n U) " “ 


- ^ sinh Oq cos p n + 0n cosh an sin P n ^ ( cos 0n£ 


3 n / 

cosh a^l) + 0 n — fl + p— sin p n sinh 0 ^ + - 

cos g n cosh a^sinh a^l - ^1 - ~ sin 3n sinh an + 

0 n cos 3n cosh a^sin f3 n £ 




The coefficients Xy n (|) and (i) are seen to be the natural mode 

shapes associated with the frequency k n . As a check on the validity of 
the Laplace transform procedure these quantities may be shown to satisfy 
the differential equations (Bl) and the boundary conditions (BIO). 

Substituting y(|,T) and from equations (Bl8) and (B19) 

into equations (B2) and (B3) gives, for the moment and shear. 


c 

M(6,t) = -2 Y_ x H n ^) sin V 

1 n=l 


(B20) 
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where 


%Ci) = 


-(i)V 


Mn 


( a n 

1 + p- sin p n sinh a n _+ 


cosh Oq cos 3 n ^ cosh a^l - ^1 - ^ sin p n sinh 
0n cosh cos p n j cos p n | - sinh cos p n + 
0 n cosh 0 ^ sin P n j^sin p n £ + ^ sinh 


and 


v (§j t ) = -2R — y X v ^(|)sin k^T 


(B21) 


where 


V 


*n 


Pn« 


n w n 


P n / Gt li 1 \ 

l^iV + Pn S±tl Pn S±Ilh 0/11 + ¥ri COS Pn C ° Sh ^) sirih “nS + 


^ ^ sin Pn slnh “n + 0n cos Pn cosh ct-nWn p n £ - 


^ sinh cos 0 n + 0 n cosh cc^ sin P n j^cos p n i + cosh cfclj 


The first eight terms in equations (£20) and (B21) have been used 
to compute the quantities M(0 ,t) and V(0 ,t). These results have been 
plotted in the range 0 < r ^ 10 in figures 4(a) and 4(b). 

Simply supported beam with an applied step end moment .- If a step 

moment M = 1 is applied to the end 5=0 of a simply supported beam 
at time t = 0, the boundary conditions are 
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kl 


?(0,t) = 0 

yU, t ) = 0 

(0 j t) = -1 
= 0 

These equations transform to (transform 7 > table I) 


(B22) 


Y(0,p) = 0 
Y(l,p) = 0 

f^(0,p) = - | 

= 0 


“'v 




(B23) 


If conditions (B23) are used to evaluate the functions Cq , C , 
Cp and in equations (B8) and (B9), the quantities Y(l,p) 

and $(£,p) become 


Y(£,p) = — -(cos mq£ - cot mq sin mq£ - 

P( m l 2 + m 2 2 ) 

cosh mg I + coth mg sinh nigl) (B24) 


$(i,p) = - 


p(mq 2 + m 2 2 ) 


“i 2 + 

— (sin mq| + cot mq cos mq| J + 


"2 2 - (^f p2 , 

^ (sinh mg I - coth mg cosh niglj 


(B25) 
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or they may be written in the forms 


Y(t,p) 


%u,p) 

Dy(p) 


(B2 6) 


where 

Ny(£^p) = sin my sinh mg^cos my£ - cosh mg£^ " cos m i si 11 ^ 1 nig s ^ n m y£ + 
sin my cosh mg sinh mgl 


and 


where 


$(^P) 


Nf(l*P) 

Ef(p) 


(B27) 


N 


$( l>p) = 


mg - 


2 - P 2 siri my (cosh mg coBh mg | - sinh mg sinh mg£^ - 


mg 


( c l\ 2 2 / \ 

Im-j^ + I — ) p sinh mg^cos my cos my|_ +- sin my sin rn^lj 


D$(p) = m 1 m 2 D y (p) 

Here again, the functions Y(l,p) and $(£,p) are single -valued, 
and the inverse transforms y(£,-r) and t(£,-r) of the quantities Yl(£,p) 
and i,p) are taken as the sums of the residues at' the poles of 

e pT Y(!,p) and eP T $( £ ,p), respectively. Consider all the roots of the 
denominators Dy(p) and Djjr(p). For my = ±mt or mg = +inn, 
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^3 


n = If 2, 3) • • • 00 (any of which, are roots of "both denominators), 
equation (B 6) becomes 


+ n 2 * 2 



+ 



(n = 1, 2, 3, . . . 00 ) (B28) 


The solutions of equation (B28) are p = -ia^ and ±ib n , 
n = 1, 2, 3; • • • 00 , where 



and these points are seen to he the locations of simple poles of both 

2/R ^ 

Y(|,p) and $(£,p). The roots p = ± -r- of the equation 

/ v O 

1 

2 2 W. 

m l + mQ =0 are zeros of the denominator, but these are also roots 


of the numerators Ny(5,p) and N^(l,p) and are not poles of either 
Y(l,p) or $(i,p). The root p = 0, on the other hand, is a pole of 
both Y(£,p) and $(g,p) for, although both numerators also go to 
zero at this point, the denominators vanish more rapidly. The equation 
m^ = 0 has three roots: p = 0, which has already been discussed,, and 


p = ±i \/R. The latter points are not poles of Y( £,p) since they are 
roots of Ny( 5,p) = 0, but they are simple poles of $(i,p). 


First consider the residues of e® T Y(t,p) and e^ T f(|,p) at the 
poles -ia n and -ib n . The sums of the residues at these poles provide 
the results 



= y~ 

n=l 

NyUjiaJ ia ^ T Ny(^-ian) 

By’(ian) Dy'(-ian) 

and 



N-^(£>ian) dan T ( ®^(£,-iaji) 

Dj'Uan) 6 Df'(-ia n r 

— 


V(ibn) 




(B29) 


D$' 


(B30) 


P" 

p- 


where the subscripts on y and i]r indicate that they constitute only parte - so-called "dynamic" 
parts - of the complete solutions for y(£,T) and ^£, 7 ). In order to obtain the values of 
Ny Ntf 

a t p = *ia n and lib^, it is necessary to use the fact that, when p equals one 

of the roots ±ia n and ilb n , one of four equations = mr, = ~nit, = inn, or = -inn 

is satisfied. The question of which value of p corresponds to each equation need not be answered 
since the following relationships can be derived: 


m jm 




mt sin nn£ 


-nit 


iPl / mo=init ‘ W'/mp =-init 2 a n 2 b n 2 + P 2 ^ 2 + b n 2 ) 


(B31) 


\ /^? \ 

fj Wy 

1 ' m-^=nji ' * ' 


-mt 


■ft) _-&) . 

' “ ' m2=inJt ' “ 'n^-inn 


(If 

n2 ” 2 - (|) 2p2 

COB 

mt£ 

2a n 2 b n 2 + p 2 | 

+ b 2 1 
+ u n j 

I 


(B32) 
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As an example, let p = ia n . Then one of the four possibilities 
= nir, m-]_ = -nir, = inn, and m 2 = -imr must occur. But for any 


of these possibilities, the ratios 


have the forms given 


and (B32) . 

Thus 

Ny(5,ia n ) 

/ C2\^ 
\%/ : 

Dy ' ( ia n ) 

•A 

^■ijr( ^ ) ia n) 

w 

Df ’(ian) 


cos nn£ 


a n 2 " b n 2 


(B33) 


The values of the other necessary ratios can be found in a similar manner 
and equations (B29) and (B30) reduce to 


o 00 

— ft \ „/ c 2 \ \ nn sin nn| /c 

yd(l ' T> =2 k) 


cos a n T cos b n T' 

a 2 V 2 


(B3*0 


td 1 (^ T ) = 


cos nit | 

7 2 _ b 2 

D n 


- 1 cos a n T 


- 1 cos b n T (B35) 


Next, consider the contributions of the pole at p = 0. Equa- 
tions (B24) and (B25) may be expanded in a Laurent series about p = 0. 
The resulting equations for Y( |,p) and f( i,p) may be written 


y(i,p) = + i + 0 (p) 

P 6 23 


t u,p) = |4 _l+ i + i + 0(p) 
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where 0(p) signifies terms of order p or higher. Thus the quantities 
e P T y(|^p) e pT f( 4,p) have simple poles at p = 0 with the residues 

^ + f (MS) 

and 

*8tU> T ) = f- - 1 + J (B39) 

respectively. It can he readily verified that equations (B 38 ) and (B39) 
actually constitute the solution of the problem when the beam is con- 
sidered to he loaded statically. 

Lastly, consider the points p = ±i\/R, which are simple poles of 
fU,p). It' can he shown that 


-M 1 

lim j, ' T = - — 
p— *-±i/R % 2R 


(B4o) 


so that the sum of the residues of e^ijKijp) atrthese poles provides 
the additional dynamic contribution 


td 2 (^ T ) = - J cos i|i T 


(B4i) 


Summing the contributions of all the poles gives 


y(S> T ) 


*(i, T ) 


= il.i! + i + 




rut sin nni /cos a n T cob b n T' 




2 -»„ s 


■( : 


= £_ 


00 

1 + i. + i-i c “' BT+8 Z 


cos nrt| 


*n 2 


(B42) 


an 2 - b n 2 |\n 2 Jt 2 


- 1 cos a n T 


Vn 2 n^ 


- 1 cos h n T 


(B43) 
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\ 

The responses y(l, T ) and i|t(!,t) are seen to have the same form as 
solutions obtained by Williams ' method (ref. 2); that is, they are the 
sum of a static part, yst(£> T ) and and a dynamic correc- 
tion, y^C S, T ) and ^a-^(?> T ) + , which is expanded in a series 

of the natural mode shapes. The y- and \|r- components of the natural 
modes have the shapes sin niti' and cos nirt, respectively, and there 
are seen to be two natural frequencies, a n and b n ‘, associated with 
each integer n so that each integer corresponds to two modes. From 
equations (B2) and (B3) , 


M(!,t) =l-|+2 


rut sin nir| 


n=l 


b n 2 


n 


i 2 * 2 


1 cos a^r - [- 5-5 - 1 cos b n Tj 


n 2 * 2 


(b44) 


and 


V( ?,t) = cos \/R t - 1 + 2R 


OO 

V cos n«| , 

Z_ a n 2 " b n 2 
n=l 


cos a^r - cos b n T) (B45) 


Six terms (n = 3) have been used in the summations of equations (B44) 
and (b 45) in computing the quantities M(i/2,t) and "V(0 ,t). These 
results have been plotted in figures 5(a) and 5(b) up to t = 8. 


Simply supported beam with an applied ramp-platform end moment . - 
Let the response to unit step moment at £ =,0, obtained above, be 
designated M]_(£,t) and V^_(I, T ). Then the response to an arbitrary 
applied moment M(0,t) may be obtained from Duhamel's superposition 
integral as 


M( I, t) 


= M(0,0)M 1 ( S,t) 


+ 



m 1 ( I, T_e)M e (o,e)d0 


v(I,t) -m(o,o)v 1 (s,t) + f V 1 (i,T-0)M e (O,e)d9 

^0 


> (B 46) 
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where the subscript 6 indicates differentiation with respect to that 
variable. Let the applied moment have the ramp-platform time history 


M(0,t) = T 
M(0,t) = 1 


Then equations (B46) become 

M(6,t) 


= J M x ( | ,T -I 


e)d e 

T-0)d0 


(0^T$ 1)_ 

( T > 1) 


(o £ t £ i) 

(t > l) 


(b4t) 


(b48) 


and 


V(8,t) 


-JN*' 


T-e)d 9 
r-9)dB 


(o £ t £ i) 

(t > 1) 


> (b49) 


Substituting and V-l from equations (b 44) and (b 45) gives for M 

and V 


(1 - 1 )t 

00 

\ nit sin nn£ 
+ p / 

/'■n 2 

n \sin a n T 

n=l a 

n 2 - *n 2 j 

\n 2 Jt 2 

/ an 

OJ 

\ sin b n T 



(0 

\n 2 n 2 

1 J~ 




H(|,t) =1-| +2 


nit sin nw| 


n=l a n 2 - b n 2 


/ b n 2 _ 
\n 2 rt 2 


(0 £ t £ 1) 


\ sin a n T - sin a n (T - l) 
/ a n 


> 


OJ 

\ sin b n T - sin b n (T - l) 

Vn 2 * 2 

7 bn J 


( T > 1) 


(B50) 
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b9 


and 


v(!,t) 





As in the preceding case, computations using terms up to_ n = 3 in 
equations (B50) and (B51) have heen made for M(i/2,t) and v(0,t) in 
the range 0 ^ r ^ 8. The results are plotted in figures 6(a) and 6(h) . 
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APPENDIX C 


EXACT CLOSED SOLUTIONS FOR UNIFORM BEAMS WITH c ± = c 2 

Introduction 


In this appendix some exact closed solutions of Timoshenko's equa- 
tions are derived for infinitely long uniform beams for which c^_ = c 2 - 
The results for beams of infinite length are then utilized to obtain 
results for beams of finite length. This is done either by simply 
restricting the infinite-beam solution to - a time interval in which it 
coincides with the corresponding finite-beam solution or by superposing 
infinite-beam solutions in such a way that the result/ satisfies the 
boundary conditions for a finite beam. 


Equations — 

If the propagation velocities c^ and c 2 are taken to be equal 
and if Tj = Oj equations (l4) reduce to 


03| + M T = 0 

(Cla) 

V| -■ 4x 2 v t = 0 

(Clb) 

0 
II 

l> 

1 

t- 

13 

+ 

jxn 

IS 

(Clc) 

ft V T - CO. = 0 

1 4x 2 

(Cld) 


where X = 


Sr, 


This system may be further reduced by writing a linear 


combination of. equations (Cla) and (Clb) and a linear combination of 
(Clc) and (Cld) to obtain 


Z| + 2X1X t = 0 


(C2a) 


Z T + 2X.iK| - 2X1Z = 0 

in terms of the complex variables Z = V +- 2X,lxo and K = M + 24iv. 


(C2b) 
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For a beam initially at rest, equations (C2) may be transformed to 
(transform 6, table I) 

+ 2 lip k = 0 (C3a) 

21i K? + (p - 21i) £ = 0 (C3b) 


in terms of the Laplace transforms £(S,p) = I e pT Z(£,T)dr and 

Jo 

poo 

k( £,p) = I e -I>T K( i,T")dT. Eliminating k produces the equation 
J 0 

£gg - p(p - 21i)£ = 0 

which has the solution 


£(|,p) = A(p)e 


-|^p(p-2X.i) 


B(p)e 


I \/p(p-21i) 


(C4) 


where A and B are governed by the boundary conditions. If the beam 
extends to infinity In the positive |-direction, boundary conditions 
stipulating that the beam remain undisturbed at | = oo must be satisfied. 
Thus, equation (C4) reduces to 

£(6,P) = A(p)e" ll/p(P " 2Xi) (C5) 

Substituting this expression into equation (C3a) produces 


*(S,p) = |^^(P - 211) e “*lf p(p - 2Xl) (C 6 ) 

The quantity A(p) may now be determined from the remaining boundary 
conditions in conjunction with equation (C5) or (C6), or both. 


Specific Problems 

Infinite cantilever given a step velocity at the root .- If the 
end 1=0 of an infinitely long beam is restrained from rotating and 
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is given a unit step velocity at time t = 0, the boundary conditions 
are v(0,r) =1 and Zd(0,t) =0. In terms of Z and K these condi- 
tions may be written 


i.p.(z(o,t£| = 0 


and 


I .P ■ (k(0,t)] = 2X 

where I.P. designates the imaginary part of the quantity in brackets. 
The boundary conditions may be transformed to (transform 7 , table i) 


and 


I.P. [1(0, pj] = 0 

i.p. g(o,p3 = ^ 


respectively. 


(C7a) 


(C7b) 


Because of the nature of these conditions, it will be convenient 
to proceed in the following manner to determine the function A(p). From 
equation (C5), it may be seen that A(p) = £(0,p); thus condition (C7a) 
establishes at once that A(p) for this problem is a real quantity. 

At |=0 equation (C6) becomes 


*(0,p) = ^|j-^p(p - 2X1) 
or, in view of condition (C7b), 


R.p.&(o,p 2 + = §XijJ p(p " 2 ^ i} (c8) 

where R.P. [^c(0,pr] is the real-part of k(0,p) . Since A(p) and p 
are real, the conjugate equation, which must also hold, i6 

R.P. f(0,pj] - 2X1 ) (C9) 



MCA TN 287^ 


53 


Eliminating R.P. [jc (0,p)J 


between equations (C8) and (C9) results in 


A(p) 



(CIO) 


If A(p) is substituted into equation (C5), the result may be 
written in the form 


e ~i\fp( P-2Xi) 

\/p(p - 2Xi) 

which has the inverse transform (transforms 1, 3> 1^ and 15 > table I) 


£U,p) = - to. 2 + 2Xi(p - <Jp 2 + bX 2 ) 


Z( 


|,t) . -4x4e iW Jo (xP^I5) - 


PT 




(r > I) (Cll) 


where J 0 and Jj_ are Bessel functions of the first kind. From equa- 
tion (Cll) the shear and angular velocity are obtained as 


V(£,t) = -kX 2 jcos Xt Jq(\Jt 2 - £ 2 ) + 


PT 


|2A(t - ©TJ 

sin xe j n l xie* - r) T — g 


T 0 (xfe 


d0| 

(t > £) (C12) 


and 


Zd(|,t) = -2.X jsin Xr J 0 (k^ 2 - I 2 ) - 


PT 


COB 


™ J 0 (^\/g 2 - ^ 2 ) Jl ^ g -' — 


dSi 


(t > 6) (C13) 
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in terms of integrals that apparently defy evaluation by analytical 
methods. 

The time history of the shear V(0,t) at the root of the beam 
may be obtained as a special case of equation (CIS); however, it will 

be noted that l{v(0,t)} = £(0,p) = A(p) (where l{v(0,t)} denotes the 

Laplace_transform of V(0,r)) so that the inverse of equation (CIO) 
gives V(0,t) directly. Thus the shear at the root is found explicitly 
in terms of tabulated functions as (transform 13, table I) 


V(0,t) = -4\ 2 [cos Xt Jq(Xt) + sin Xt Ji(\t)T] (r > 0) (Cl4) 

Since A(p) may be replaced by l{v(0,t)} , the inverse transform 
of equation (C6) may be written for this problem in the form (trans- 
forms 1, 2, 3, and 17, table I) 


r 

''T 

1 e ixl v(o,t-8) - 

l 

V(O,T-0)e i>1 ^ 

£ 


es 1 




- £ 2 


ij 0 UJe 2 - I 2 ) 


d0 


(x > |) 


or, after integration by parts., 
K(£,t) = 2Xie j 


i iXT J 0 {x\It 2 - | 2 ) - jv(O,T-0) + 

T 2 ^ i Q . sin X(t - 0 ) J x [x(t - e3je 1X ^Jo(x/e 2 - £ 2 ) 


d0 

(t > I) (C15) 


where V(0 ,t) is defined by equation (Cl4). From equation (C15) the 
moment and linear velocity may be written as 


M(|,t) = -2.X sin Xt Jo(x^t 2 - i 2 1 - 2.X 


> - » rft 


V(o,T-e)cos x© - 


^ f sin xfy - e) j x [x(t - J 0 (^/e 2 - (t > l) (ci6) 
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and 


v(£,r) = cos \t J 0 (x/r 2 - i 2 ) - f— V(0,T-e)sin 10 + 

^ 1 2A, 

~~^ G - sin l(r - e) J X [X(T - Q)^J 0 {x\1q 2 -l 2 )d 9 (r > |) (C17) 


again in terms of integrals requiring evaluation by approximate methods. 

As in the ca_se of the shear, it will be more convenient to deter- 
mine the moment M(0 ,t) at the root by means other than a reduction of 
of the general relation, equation (Cl6) . Since, for this problem, 

R.P. [k(0,p)] = l{m(0,t)}, equation (c8) may be written, after substitu- 
tion for A(p), in the form 


l{m(0,t)} = M'Jp 2 + bX 2 - p) 


The inverse transform of this equation is (transforms and l4, table i) 

pr 


M(0,t) = 21 


Jl(2ie) 


0 


cL0 


(t > 0) (Cl8) 


But, from reference 17, 


^ = 2XJ (21t) - ^[ji(2XT)] 


and 


J o (2ie)d0 = t|j 0 (21t) + J[j 1 (21t)H 0 (21t) - J 0 (21t)h 1 (21t)]| 

where Hq and are Struve functions. Thus, equation (Cl8) becomes 

M(0,t) = 21-|x.t[ 2 - nHi(2\TT] J 0 (21t) - [l - t[1tH 0 (21t)] J^(21t)|» 


(t > 0) (C19) 
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The Struve functions Hq and H]_ are tabulated in reference 17 
in the range 0 ^ 2XT ^ 15 *9* For larger values of the argument, it is 
convenient to use the approximations (ref. 17) 

H 0 ( 2 Xt) ss Y q ( 2 Xt) + JL ( 2 Xr > 15 . 9 ) 

H 1 (2Xt) ^ Y 1 (2Xt) + i^l + 

where Yq and Y^ are Bessel functions of the second kind. From 
reference l 8 , 

Y 0 (2Xr)J 1 (2Xr) - Y 1 ( 2 Xt)J q ( 2 Xt) = JL 

so that 

M(0,t) 2X [l - -° - gxV T) J (2 XT > 15.9) (C20) 

Computed values of V(0,t) and M(0,t) , obtained from equa- 
tions (Cl4), (C19) , and (C20) for a beam with X =• 5? have been plotted 
in figures 4(a) and 4(b), respectively, in the range oSt =2. In 
this range, the root of a finite cantilever beam behaves as if the beam 
were infinite, since the effect of the free end is not felt until the 
return of the wave front at t = 2 . 

Infinite simply supported beam with an applied step end moment . - 
If the end £ = 0 of an infinitely long beam is simply supported and 
is subjected to a unit step bending moment at time t = 0 , the boundary 
conditions are M(0 ,t) = 1 and v(0,t) = 0, 'or K(0 ,t) = 1. This 
condition transforms to (transform 7 , table I) 

•s(0,p) = i ■ (C21) 

Substituting equation (C 6 ) into condition (C2l) reveals that 


kX e ) 


(2Xt > 15-9) 


A(p) 


2Xi 

\/p(p - 2 Xi) 


(C22) 
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so that equations (C5) and (c6) become 


£( S,p) 


e -l\/p(p-2X.i) 
\/p(p - 2X1) 


(C23) 


and 


K(S,P) = 1 e- S ^ (p - 2Xi> (024) 

P 

respectively. The inverse transform of equation (C23) is (transforms 3 
and 15 , table i) 


Z(£,t) = 2Xie 1X,T J 0 (x^r 2 - I 2 ) (r > £) (C25) 

and the inverse of equation (C24) may be written (transforms 2, 3, **■, 
and 1 6 , table I) as 


K(|,r) = e ±Xi 





1X8 

e 


ji {^ z - 1 2 ) 
W - s 2 


d 0 


(t > |) (C2 6 ) 


From these equations, the shear, angular velocity, moment, and linear 
velocity are obtained as 


V(|,t) = -2X. sin Xt J 0 [x\ir 2 - £ 2 ) 

co(£,t) = cos Xt jJ^Xyfr^ - | 2 ) 

Jl(x/e 2 - £ 2 ) 


T 

M({-,t) = cos X.| - X£ cos xe i - 

v/o2 _ e c: 


v(£,t) = sin - X| J sin X6 


/© - i 

Ji(4 


e 2 - § 2 ) 


'e 2 - l 2 


d0 


d0 


(t > I) (C27) 
(t > £) (C28) 
(t > |) (C29) 

(t > 6) (C30) 


Thus, V" and 55 are determined everywhere in closed form in terms of 
tabulated functions, but the relations for M and v contain integrals 
which must be evaluated by approximate methods. 
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The results of this section are utilized in the next where a finite 
simply supported "beam is considered. 


Finite simply supported “beam with an applied step end moment .- If 
the end 1 = 0 of a finite simply supported beam is subjected to a 
unit step bending moment at- time t = 0, the boundary conditions may be 
written 


M(0,t) = 1 

v(0,t) = M(1,t) = v(1,t) = 0 


or, in terms of K, 


K(0,t) = 1 
K(1,t) = 0 


(C31) 


(C32) 


Although a direct solution of equations (C2) is possible in this case, 
the response of the finite beam may be obtained by a somewhat simpler 
procedure in which the responses of infinite beams to the same disturb- 
ance are superposed. This procedure is described below. 


A series of semi-infinite simply supported uniform beams extending 
in opposite directions are shown in figure 7* The beams have been posi- 
tioned in the figure so that the origins of the space coordinate I lie 
on the same vertical line, and the segments of greatest interest, 

0 | ^ 1, have been outlined with solid lines while the rest— of each 

beam is defined by dashed lines. For each beam the origin of the 
coordinate | is seen to lie at a different position relative to the 
end where a unit positive moment is suddenly applied at time t = 0. 

An infinite number of these beams Is assumed to exist. 


The response of the top beam has already been determined as 
(eqs. (C25) and (C26)) 


Z(l,r) 


2XIe 1Xx J 0 (\y/T 2 - § 2 ) 


(t > I) 


and 


K( I, t) = G( |, t) 


(t > 6) 
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where 


G( t, T ) 



The response of the second beam may be determined from these equations 
simply by accounting for the changed direction and shifted origin of 
the coordinate i . This process gives 


Z(i/r) = 21ie iXr J 0 [x^T 2 - (2 - |) 2 J (t > 2 - I) 

and 

K(6,t) = -G(2-J,t) (t > 2 - |) 

In order to obtain the response of the third beam from the response of 
the first, only the shifted origin need be accounted for. Thus, for 
this beam. 


Z(S,t) 


2\ie 1XT J 0 X\It 2 - (2 + |) 2 


(t > 2 + i) 


K(S,t) = G(2+|,t) (t > 2 + |) 

Similarly, for the fourth beam, 

Z(S,t) = 21ie 1XT J 0 [ljT 2 - (4 - j (t > k - 6) 

K(|,t) = -G(4-|,t) (t > 4 - I) 

and so forth, for all the rest. In each case, the region where the 
response is not specified is a region of zero response. 

Let the responses of all these beams to their respective dis- 
turbances be superposed and consider only the response of the segment 
0 ^ i ^ 1 of the resulting composite beam. Since the wave fronts of 

all the disturbances travel at the same velocity ^ = 1, the response 

of this segment may be written 
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Z(£,t) 

Z(S,t) 

z(6,t) 

Z(|,t) 

and 



(C33) 


K(|,t) = 0 
K(£,t) = G(|,t) 

K(£,t) = 0(£,t) - G(2-6,t) 

K(5,t) = G(|,t) - G(2-6,t) + G( 2+|, t) 


(O £ T H £) 

(I ^ T ^ 2 - £) 

(2 - i ^ T < 2 + I) r (C34) 
(2+|^t^4-£) 


Since G(0 , t) = X, equations (C34) are seen to satisfy the boundary 
conditions (C32) and the response obtained above must be that of a 
simply supported uniform beam to a unit dimensionless step moment applied 
at the support | = 0 . 
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The vertical shear V and bending moment M may he obtained from 
equations (C33) and (C3^) as 


V(S,T) = 0 

V(i,T) = -2\ sin XT Jq(x\It 2 - | 2 ) 

f 

V(£,t) = -2.X sin Xr 


(0 < r< £) 
(I g T < 2 - I) 


|j 0 (x.f 2 - * 


J Q \x{i 

|jo(i'/ t2 - 1 2 )+ j 0 E ^ 2 - ( 2 ■ + 

j 0 [m/t 2 - (2 + 


J 0 \x\jr 2 - (2 - t) 
V(£,t) = -2X sin XT 


(2 - I ^ T £ 2 + 6) >(C35) 


(2 + £^t^4 - 6) 




and 


M( I, t) =0 
M(6,t) - H(|,r) 

M( S,t) = H( 6,t) - 1(2-1, t) 


(0 ^ |) 


(i = t = 2 ~ l) 


(2 - I < T < 2 + I) MC3 6 ) 


= H( I, t) - H(2-5,t) + H(2+I,T) 


(2+ I^t^4 - I) 
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where 




= cos X| - X| 



The angular and linear velocities oj and v may, of course, also he 
obtained from equations (C33) and (C34) . 

Computations have been made for the quantity V(0 ,t) from equa- 
tion (C35) with X = 5* The results have been plotted in figure 5(a) 
in the range 0 ^ r ^8. 

Finite simply supported beam with an applied ramp-platform end 
moment .- Let the applied moment at the end 1 = 0 of a simply sup- 
ported uniform beam have the time history 


M(0,t) = T 

(0 < T < 1) 

H(0,t) = 1 

(t > 1) 


T 




(C37) 


The response of the beam to this disturbance may be obtained from the 
response to a unit step by using Duhamel's superposition integral. 
Thus, if the response to a unit step, as given by equation (C33)> is 
designated Zi(|,t), the response Z(£,t) to the applied moment (C37) 
may be written 


Z(S,T) 



ZiCi^-ejde 


(0 ^ T ^ i)"l 


z(I,t) 


z 1 (i,T-e)ae 


> (C 38 ) 


(T " X) J 
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Substituting from equations (C 33 ) and letting i = 0 results in the 
following expressions for the end response Z( 0 ,t): 


Z( 0 ,t) = 2 \i 


r 

'0 


e iX(T-e)j o j^( T _ ©yj a.© 


Z( 0 ,T) = 2 x 1 J' e 1 X ( T-e )j 0 [x.(T - ©Tj de 


Z( 


Z( 


Ij-Xi e 1 ^ T_e ^J 0 |x 1 /(T - Q ) 2 - k 
0 , t ) = 2 Xi f 1 e^< T -«)j 0 [l(T . 8 )| d© + 


4Xi e iA.(T-©) J( , 

Jo 


x^(t - 9 ) 2 _ i(. ae + 


(0 £ T £ l) 


( 1 = T = 2 ) 


0,t) = 2Xi J' e lx ( T " e )j 0 |y,(T - ©7J d© + 

4Xi e iX(T -fl)j 0 j\^( T - ©) 2 - ^jd© (2 <t^3) 


Z( 0 ,t) = 2 Xi J' e 1 X ( T " 0 )j o |x(T - © 7 ]d© + 


/■( c 39 ) 


(3 ^ r < k) 


iai e 1 X ( T ~ & )j 0 X\j ( t - ©) 2 - 16 d© ( ^- = t = 5 ) 


J 



6b 


NACA TN 2874 


It can be seen that (transforms 3> b, and 9, table i) 


L 1 


e iX0 j o (\e)dei> = 1 


i d / I P - 2Xi 


P y P (p - 2Xi) X- d P \ V P 
or (transforms 5 and 11, table I) 


- 1 


-'T 

Jo 


e 1Xe J o (X9)d0 = re iXT [j 0 (Xt) - iJ-^X-r^ 


Then, Z(0 ,t) may be reduced In the region 0 ^ r ^ 2 to 


Z(0,t) = 2X1 jT e iX0 j o (xe)de 


Z( 0 ,t) = 2XiTe iXr [ j 0 (Xt) - iJ^X-rj] (0 g t g l) (c40) 


and 


Z(0,t) = 2Xi e 1X0 J o (X6)dfi 


= 2Xi 


e iXe «r o (Xe)d0 - e iX0 J o (X6)d0 


Z(0,t) = 2XiTe lXT [j 0 (XT) -IJ^Xt)] - 

2Xi(T - i) e iX ( T - 1 )|j 0 [x(T - l)\ - iJifxfT - 1)]J 

(1 ^ t % 2) ( Cbl ) 
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From these equations, the shear at 1=0 may he written in this range 
as 

A 

?(0,t) = 2Xt[cos Xt ^(Xt) - sin Xr J 0 (Xt)] (0 ^ t £ l) 

V(0,t) = 2Xt£cos Xt J]_(Xt) - sin Xt Jq(Xt)J - 


2X(t - l)|cos X(r - D^pr - l)] - 
sin X(t - l) Jq [x(t - l)]^- 


> (Cl+2) 


(l ^ T 


— 2) 


J 


In the range t >2, however, the integrals in equations (C39) apparently 
cannot he evaluated analytically. If the shear Vp(0,T) due to a unit 
step moment has heen computed, it may he convenient to write, for the 
shear at the point £ = 0 resulting from the disturbance (C37) , 


V(0,r) 



V 1 (O,T-0)d0 


V(0,t) 


~1 _ 

V-lCO/t-S)^ 

Jo 


(0 ■£ T £ 1 ) 


(t > 1 ) 


> (C43) 


and then evaluate the integral V]_(0 ,t -6) d6 numerically to obtain 

V(0,t) in the range r > 2. 

The quantity V(0 ,t), for a beam with X = 5> has heen computed in 
the range 0 ^ t % 2 from equations (C^2) and has heen obtained in the 
range 2 ^ t % 8‘ by a numerical integration of the exact curve of 
figure 5(a) in accordance with equation (0^3). These results are 
presented in figure 6(a) . 
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TABLE I.- LAPLACE TRANSFORMS 
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(a) Space -time plane. 

Figure 3.- Characteristic grid for application of numerical procedure to 
uniform beams with disturbance applied at 5 = 0. c-^ = Cq- 


(a) Time history of the nondimensional vertical shear at the root. 

Figure 4.- Response of a uniform cantilever beam subjected to a step 

velocity at the root. 








Period of first 
natural mode 






(b) Time history of the nondimens ional moment at the center. 

Figure 5.- Concluded. 
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Time, r 


(a) Time history of the nondimens ional vertical shear 

at the end | = 0 . 


Moment, M, 
at 



(b) Time history of the nondimens ional moment at the center. 

Figure 6.- Response of p.'. uniform simply supported beam subjected to a 

ramp -platform moment at | = 0. 
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Figure 7- - Semi-infinite uniform beams su 
of a finite simply supp 
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